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Abstract 

Given a bidirected ring with capacities and a demand graph, we present an approxima- 
tion algorithm to the problem of finding the minimum a such that there exists a feasible 
unsplittable routing of the demands after multiplying each capacity by a. We also give an 
approximation scheme to the problem. 

1 Introduction 

In this paper a bidirected ring means the union of two oppositely directed circuits on the same 
set of nodes. The motivation of the present paper is to compute a routing of demands in an SDH 
(Synchronous Digital Hierarchy) based network. Such a network is a transmission network and 
we will refer to it as the supply graph. Also, we will refer to the set of demands as the demand 
graph. The backbone networks of the European telecommunication companies are based on 
the SDH technology. The ring architecture plays a key role in this technology, because it 
provides an efficient self-healing mechanism in case of failure: although the rings are bidirected, 
only one direction is used, while the other would be useful in case of a failure. Most of the 
SDH backbone networks consist of such uni-directed circuits. From such uni-directed circuits 
different topologies have been built, like tree of circuits or cycle of circuits. In the next section 
we will show that the routing in a supply graph which is a cycle of circuits can be reduced to a 
routing problem in a bidirected ring. One can find more about the SDH transmission protocol 

Formally, we are given a directed supply graph G with capacity function c: E{G) M+ 
and a directed demand graph H on the same node set with demand values d: E{H) R+. A 
routing of is a collection of Mf-paths of value fc, one for each demand edge of E{H) joining 
u io V with value k, satisfying the capacity constraint c. Since each demand is routed along 
a unique path we will speak about unsplittable routings. In this paper a routing is always 
unsplittable, unless when speaking about the fractional solution of the linear relaxation of the 
problem. We say that a routing of H has load at most (less than) c if for each edge e € E{G) 
the sum of the values of the paths using e is at most (less than) c(e). 

In this paper we consider the routing problem in bidirected rings. Given a bidirected ring G 
with capacity function c : E{G) — >■ M+, and a directed demand graph H with demand values 
d: E(H) K_|_, we want to find the minimum a for which there exists a routing of H with 
load at most ac. This is called the balanced bidirected ring routing problem. Under 
balancedness we mean that the value of a is introduced to ensure that the remaining capacities 
are as large as possible, allowing the network to carry larger demands. The case when c is 
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uniform (that is, all capacities are equal) and each demand is 1 was solved with an elegant 
method by Wilfong and Winkler 15 . In Section |3] we follow the lines of their proof to give an 
algorithm which returns a routing with load less than aoptc + |-D where aopt is the optimum 
solution of the problem and D is the maximum value of the demands. We also show that this 
error term can indeed occur. 

In Section 31 using a method of Khanna [2], we show an approximation scheme to the 
problem which for any e > yields a routing with load less than aopt{c + e^^), with a trade-off 
with the running time. Here n — \V{G)\ and X)"^ t^*^ the capacities of the edges. 

By solving a linear program, a fractional solution of the balanced bidirected ring 
ROUTING PROBLEM Can be given, that is, we can calculate the minimum a, denoted by a*, 
such that there exists a fractional routing of H with load at most ac. Clearly a* < aopt- 
Our approximation algorithm for the balanced bidirected ring routing problem gives 
a routing such that its load is actually less than a*c + If a* < 1 then our solution may be 
considered good in practice since it requires less than |D additional capacity on each edge. On 
the other hand, if a* > 1 then there exists no feasible routing with the given capacities. Hence 
we are facing a network design problem, and our goal is to increase the capacity function with 
minimum cost, in order to have a routing satisfying the increased capacity. We show in Section 
Elthat even this problem can be handled with techniques similar to Section |3l 

If the capacity function c is uniform and the ring is undirected, the balanced ring routing 
problem was studied by many researchers due to its significance in telecommunication net- 
works. This problem is called the undirected ring loading problem, first considered by 
Cosares and Saniee [T]. Schrijver, Seymour, and Winkler [3] gave a combinatorial approxima- 
tion algorithm for the undirected case, if c is not necessarily uniform. Their algorithm returns 
a routing requiring less than |Z) more capacities on each edge than in an optimum solution. 
The balanced bidirected ring routing problem was first considered by Wilfong and 
Winkler [5] who gave an exact algorithm for finding an optimum routing in a bidirected ring, 
in case c is uniform and each demand is 1. Our considerations in Section [3] are based on their 
method, yielding a generalization of their result. Our result are more general than theirs in the 
aspect that in our case the demands are not restricted to be 1, and the capacity function is not 
necessarily uniform. 

We point out that the balanced bidirected ring routing problem is NP-complete. 
Indeed, the partition problem can be reduced to it in a straightforward way, just as in 
the undirected case (Cosares, Saniee [1]). Moreover, contrary to the undirected case, the cut 
condition is not sufficient for the existence of a fractional solution. That may be the reason that 
no combinatorial algorithm is known finding an optimum fractional solution of the bidirected 
ring routing problem, unlike in the undirected case. 

2 Motivation 

In some SDH backbone networks the supply graph G forms a cycle of directed circuits, with 
capacity function c : E{G) — > R+. In this section we show that the routing in G can be reduced 
to a routing problem in a bidirected ring. Two neighboring circuits of G meet in a node called 
a common node. The common nodes divide each directed circuit into two edge sets, called 
arcs. Sec Figured) 

Let / G E{H) be a demand edge with value d joining vi to V2- Let ai be the common node 
a of G minimizing distcivi, a). Similarly, let 02 be the common node a minimizing distQ^a, W2). 
If V2 is contained in the fiOi-path of G then only one ?;i?;2-path exists in G so simply delete / 
from H and decrease the capacity of this path by d. Otherwise all t;iW2-paths use the supply 
edges on the wioi- and on the a2W2-paths of G. So replace / by a new demand edge /' joining 
oi to 02 with value d and decrease the capacities of the supply edges on the viUi and on the 
a2W2-paths by d. After doing so for each demand edge, the new demand graph H' has a routing 
in G with load at most the new capacity function c' if and only if H has a routing with load at 
most c. Moreover, the non-common nodes of G are incident to no edge of H' , so we can think 
of an arc A as only one edge with capacity the minimum of c'(e) taken over all edges e E A. 
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Figure 1: Reduction to a bidirected ring 



The new supply graph is a bidirected ring with n nodes and 2n edges, where n is the number 
of circuits in G. 

3 An approximation algorithm 

In this section we show an algorithm yielding a routing in a bidirected ring with load less than 
aoptc+ |D where aopt is the solution of the balanced bidirected ring routing problem 
and D is the maximum value of the demands. Our algorithm is a modification of that of Wilfong 
and Winkler ^ who solved the balanced bidirected ring routing problem in the case 
if c is uniform and the demands are 1. We do not count running times since solving a linear 
program is included. 

From the two directions of the ring we say that one is the forward and the other one is the 
backward direction. Accordingly, an edge e 6 E{G) can be forward or backward, and from 
the two possible uw-paths {u,v G V{G)) one is the forward and the other one is the backward 
path. For an edge / £ E{H) joining u to v, an /-path refers to any of the two uu-paths of G. 
The edge sets of these two paths will be denoted by F{f) and B{f), respectively. 

The first step is to solve the LP-relaxation of the problem. There are two possibilities of 
routing demand / G E{H) hence we introduce a variable < </9(/) < 1 with the meaning that 
^p{f) fraction of the demand / is routed forward and 1 — ^p{f) fraction is routed backward. So 
we have the following LP-relaxation, whose optimum solution is denoted by a* . Both sums run 
on demand edges / G E{H). 

min a, s.t. 

0<^<1 (1) 
E vU)d{f)+ J2 (l-¥'(/))rf(/)<ac(e) WeeE{G). 

Note that we would get aopt if ^fi was required to be integer. Now we manipulate the 
demands / G E{H) with < ^p{f) < 1. Such demands are called split. We say that demands 
/ii /2 G E{H), where joins Si to ti for i = 1,2, are parallel if the end nodes are placed in the 
ring in the order si, ti, t2, S2 (some of these nodes may coincide). Assume that /i, /2 G E{H) 
are parallel split demands. Call the /^-path containing both s^^i and ^3.^ the long /^-path 
for i = 1,2. Let Xi denote the amount of flow of fi sent along the long /^-path in our fractional 
solution. If we reroute min(xi,a;2) amount of flow from the long /^-path to the other /^-path 
for i — 1,2 then one of the demands /i, /2 will not be split any more, moreover, we do not 
increase the load of any edge of G (we may even decrease it somewhere). So at most |£'(iJ)| 
such uncrossing steps are possible and flnally we get a fractional solution where there are no 
pair of parallel split demands. Especially, demands with the same source node are parallel. So 
it will hold that for each s G V{G) there exists at most one split demand / G E{H) with source 
s. 
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Figure 2: An optimum fractional routing 

Denote the nodes of G by si, . . . , s„ in the forward order. Now we try to unspht the 
remaining spht demands. Let fi be the spht demand with source Si (if any). Assume that Xi 
fraction of fi is routed forward and yi fraction backward. Let Wi = yi if we would set ip{fi) 
to 1 and Wi — —xt if would be set to 0. If we round ^pifi) to or 1 then the load of an 

edge e € E{G) increases by Wi if e is contained in the forward /i-path, it decreases by Wi if e is 
contained in the backward /^-path, and it does not change elsewhere. There are no two parallel 
split demands so the change of the load of an edge is ±X]j<i<fc ^« ^'^^ some j < k where k 
may be greater than n but then the indices of w are meant modulo n. Here the sign depends 
on whether e is a forward or a backward edge. Now we try to set the value of Wi to yi or —Xi 
for all i in such a way that | X]j<i<fc < holds for all j < k. To achieve this it is clearly 
sufficient that 

- D/2 < ^ < D/2 holds for all k < n. (2) 

l<i<k 

^ can be easily achieved by greedily setting the values Wi to yi or —Xi one after another, since 
Xi + yi < D holds. Hence by this procedure each load is increased by less than ^D. 

So the algorithm for the balanced bidirected ring routing problem is the following. 

1. Solve the LP-relaxation of the balanced bidirected ring routing problem. 

2. Uncross the parallel split demands. 

3. Unsplit the remaining split demands in the above described greedy way. 

Theorem 3.1. The algorithm gives a routing of H with load less than a*c+ < aoptC+ |-D. 

Proof. After the uncrossing procedure the routing still has load at most a*c. Moreover, as we 
observed, the load of each edge e £ E{G) increases by less than during unsplitting. □ 

Our goal was to find a routing with load at most aoptc + |-D, but actually, similarly to the 
undirected case [3], the routing we got has load less than a*c+ |D < aoptC+ |-D. However, this 
fact does not always induce improved efhciency. Indeed, next we show that for every S > there 
is an example where the output routing of the algorithm has load not less than aoptC+ — S. 
Such tight examples are very important in understanding how algorithms can be improved in 
order to decrease running time. 

Consider the capacities of the ring as shown in Figure [2] The edges without capacity, as well 
as the dashed edges have large capacities in both directions and e is a number which is small 
enough with respect to k. All demands are of value 1. Observe that the cut {e', e"} is tight, 
that is, the total value of the demands crossing it equals c(e') + c(e"). It is easy to see that 
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there is no integer routing with a = 1 due to the cut {e', e"} with fractional capacities. The 
best integer solution is obtained if the capacity of e" is raised to fc + 1, that is a^pj = ^.^t^f.^ ■ 
Indeed, e is small enough so we can route the demand of node v forward while all other demands 
backward. Thus limfe_>oo ct-opt = 1 by an appropriate choice of e = e^*^^. 

Clearly a* — \, and an optimum tp (which is a vertex of the polyhedron (U}) is shown in 
brackets in Figure[2j Assume that the unsplitting procedure starts at node u. The split demand 
at node u is routed forward. Now the split demand at node v is routed backward and then the 
demand of node w forward. All remaining demands will be routed backward, according to our 
greedy heuristic. Now e' has load 2 and so the error at e' tends to 



2 - (1 + fee) ■ a^^^ _ 3 



lim 



by choosing ke 0, and using that D = 1. Note that the fractional capacities can be made 
integer by scaling. 



4 An approximation scheme 

If the value of the demands are small with respect to the minimum load of an optimum routing 
then the |Z3 additive error is a small deviation from the optimum solution. Otherwise this error 
term can be significant. In this section we develop an approximation scheme for the balanced 
BIDIRECTED RING ROUTING PROBLEM, which for any e > gives an algorithm polynomial in 
the number of nodes, finding a routing with load less than aopt(c + ec) where n — \V{G)\^ and 

- _ — s^^^p . We use the solution method of Khanna [2] who presented an e-approximation 

scheme for the case of an undirected ring with c uniform. By definition, an e-approximation 
scheme for a minimization problem is an algorithm which, for any e > 0, returns a solution 
with value at most 1 + e times the optimum value. In this sense the method of this section 
is not an e-approximation scheme since it does not approximate aopt- However, for M > 0, 
in the class of balanced BIDIRECTED ring routing problem instances, where c < Mc(e) 
holds for all edge e, our method is an Me-approximation scheme. Many problem instances in 
practice arising from an SDH backbone network belong to such a class with reasonably small 
M, in particular the instances where c is uniform. 

For s,t e y{G) the long si-path is defined to be the longest st-path in G, or to be the 
forward one if the two paths are of equal length. Recall that a* denotes the fractional optimum 
to the balanced bidirected ring routing problem. Let a' > a* be any real number and 
let 

E'^{fe E{H) : d{f) > fea'c} . 

Note that 

Hence \E'\ < ^n. If a' > aopt then less than | edges of E' are routed in the long path in any 
optimum routing, since otherwise the sum of the loads would be more than 

lea'c ■ f • f = na'c = a' J2eeE{G) ^(e) > aopt J2eeE{G) c(e). 

which is impossible. Thus, independently of whether a' > aopt or not, for all subsets E" C E' 
with \E"\ < I we do the following. We route the demands of E" in the long paths and the 
demands of E' — E" in the short paths. Denote the load of e G E{G) in this routing of E' by 
^(e). Now denote by a*^,, the optimum of the following linear program, where the sums run on 
/ e E{H) - E'. 

min a, s.t. 
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< (/? < 1 

E E (l~^(/))</) + Ke)<ac(e) Ve G i?(G). 

e6F(/) eeB(/) 

Note that the maximum value of a demand in E{H) ~ E' is at most ^ea'c. Hence exactly as 
in the previous section, we can find a routing of E{H) — E' with load less than a*^,,c — I + ea'c. 
In the case when a' > aopt, in any optimum routing of H less than - edges of E' are routed in 
the long path, hence we get that one of the above routings has load less than aoptc + ea'c. For 

any a' , from these X]i=o ■ routings of E{H) choose the one with load at most ac + ea'c, 

such that this a is minimum. 

The number of subsets E" C E' to consider is Yll^^o {^^- • Using that < (efc)'/^' 

holds for any integers k > I, when \E'\ > - we get 




Since the left hand side is monotone increasing in \E'\, the above bound is valid when \E'\ < |, 

too. Hence the number of subsets E" C E' to try is Oe{ni). 

Now we show how a' can be chosen. It is clear that a* < aopt < '2a* , since there are exactly 
two paths between any two nodes in G. So first determine the value of a* and then run the 
above algorithm with a' — ai — -^^^a* for i = 0, . . . , for some integer N. As we mentioned, 
at the point when aopt l£ ctj !i (1 + ^/N)aopt happens to hold, our routing has load less than 
c^optic + e{l + 1/N)c). Thus, finally, from these + 1 routings choose the one with load at 
most a(c + e{l + l/N)c), such that this a is minimum. In the beginning we should also replace 
e by e/ (1 + in the above scheme. 



5 Network design 

The approximation algorithm presented in Section [3] gives a routing in a bidirectcd ring with 
load less than a*c+ Recall that we can calculate the fractional optimum a* by solving a 
linear program. Observe that a* < 1 means that there exists a fractional routing in the supply 
network G with the given capacities. Since our approximate integer solution requires less than 
|Z? additional capacity on each edge (hopefully, much less on average), this solution may be 
considered good in practice. On the other hand, if a* > 1 then there exists no feasible routing 
with the given capacities. Hence we are facing a network design problem, where we want to 
increase the capacity of some edges of the ring, with minimum cost, in order to guarantee the 
existence of a routing satisfying the increased capacities. 

Formally, a widening cost We is given on each edge e G E{G), measuring the cost of 
increasing the capacity of edge e by one unit. We aim to find the minimum of X]ee_E(G) 7e^^e 
such that there exists an integer routing with load at most ("fe + c(e))(l — a) on any edge 
e G E{G). Here a is some robustness factor known a priori. Next we give a heuristic algorithm 
for this problem, with similar methods than that of Section |3l First we solve the fractional 
relaxation of the problem. 

min E] leWe, s.t. 

eeE(G) 

< < 1 (3) 
E V{f)d{f)+ E (l-^(/))d(/)<(7e + c(e))(l-a) Ve G £;(G). 

eeF(/) e6-B(/) 
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Next uncross the split demands, and then unspht the remaining spht demands, just as 
in Section [3l Exactly as in Section |3l we get an integer routing of H with load less than 
(7* + c(e))(l - a) + ff, with the optimum 7*. 

Hence we have a routing with widening cost which is at most ^ X]eG-B(G) more expensive 
than the optimum cost. We can also choose a to ensure that the load of this routing is at most 
7* + c(e), yielding a feasible unsplit routing in the new, increased network. 

We mention that it is also possible to apply the method of Section|4]to yield for any e > an 
algorithm polynomial in n, finding a routing with widening cost at most X]eG_E(G) '"'e more 
expensive than the optimum cost, where C* is the minimum sum of loads of fractional routings 

of H in the uncapacitated ring G (one should simply define -E' = |/ G E{H) : d{f) > , 

where C* < C < 2C*). 
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